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Abstract. A synchronization procedure of the generalized type in the sense 
of Rulkov et al [Phys. Rev. E 51, 980 (1995)] is used to impose a nonlinear 
Malasoma chaotic motion on the Frenet-Serret system of vectors in the differential 
geometry of space curves. This could have applications to the mesoscopic motion 
of biological filaments. 
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The Frenet-Serret (FS) vector set is in much use whenever one focuses on 
the kinematical properties of space curves. The evolution in time of the FS 
triad is one of the most used descriptions of the motion of tubular struc- 
tures such as stiff (hard to bend) polymers pQ. Many biological polymers 
including the DNA helical molecule are stiff and their movement is of funda- 
mental interest. At the mesoscopic level many sources of noise and chaotic 
behaviour affect in a substantial way the motion of the biological polymers. 
In general, one can think that a synchronization between the motion of the 
polymers and the chaotic (or noisy) sources could be achieved in a natural 
way through some control signal. We illustrate this idea employing a gen- 
eralized synchronization procedure based on the theory of nonlinear control 
by which we generate a chaotic dynamics of the FS evolution equations 

f = k ■ N 

N = T-B-K-f 

B = —t ■ N . 



With this goal in mind, one should first write the FS system in the form 

x = Ax + EU y = Cx, (1) 

where x € B? is the vector having as components the tangent T, normal 
N and the binormal B, whereas A is the transfer matrix, E is an initial 
vector that determines the channel where the control signal is applied, and 
C = (01,02,03) determines the measured signal of the FS system. 

In this way, the main objective, from the synchronization standpoint, is 
to force the states of the slave system, which is the FS system, to follow 
the trajectories of the master system that, in general, presents a chaotic 
behaviour. This is achieved by applying, in a well-defined way, a signal 
given by U = $(x) and rewriting the FS system in the form given in Eq. Q 
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y = (ci c 2 c 3 ) N 
V B , 

where, ci = 1, c 2 = 0, C3 = 1. If we choose now B = (0 0) T as the 
input vector, one gets the open loop FR dynamics. This type of dynamics 
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is shown in Fig. 1 for k = 1, r = 0.9, and initial conditions given by To = 
0.0024, N = 0.0026, and B = 0.0039, where the states of the system are 
given by the derivatives of the tangent, normal, and binormal unit vectors, 
respectively, whereas k and r are the curvature and torsion scalar invariants. 

On the other hand, a chaotic oscilator is a dynamic system whose evolu- 
tion is difficult to predict. In general, its main feature is the sensibility to 
the initial conditions and the variations of the parameters. Thus, its long- 
term behaviour is hard to estimate. In the following, we will use one of the 
simplest chaotic oscillator systems that has been introduced by Malasoma 

m 

X\ = X 2 

X 2 = X 3 (3) 

X 3 = -01X3 — X\ + XiX 2 , 

where a is the bifurcation parameter. In matrix form, we get 
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y = (10 0).. (4) 

This system exhibits chaotic behaviour for 2.0168 < a < 2.0577. The chaotic 
evolution is shown in Fig. 2 for a = 2.025 and initial conditions X\ = 0.0022, 
X 2 = 0.0024, X 3 = 0.0039. Having the two systems in the matrix form, we 
choose the Malasoma one as the master system and the FS system as the 
slave one, i.e., the Malasoma dynamics will be imposed to the FS motion 
through the signal U = That means that a nonlinear dynamics is 

forced upon the FS system leading to its chaotic behaviour. 

To get the chaotic FS system one should achieve the synchronization 
between the master (of subindex M in the following) and the slave systems 
(of subindex S in the following) . For this, one defines a third system, which 
refers to the synchronization error given by the difference in the dynamics 
of the two systems, i.e., 

e'i = Ke 2 + Xi(X M ) 

e'2 = re 3 - Kd + X 2 (Xm) (5) 
e'3 = -re 2 + A 3 (X M ) - U , 



where 



Ai = X 2 m(1-k) 

X 2 = X 3 (1-t) + kX 1m (6) 
A3 = —aX 3 M + Xim(X 2 m — 1) + t X 2 m ■ 
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The function U = $(e) gives the control action that leads to the synchro- 
nization of the two systems. 

Once (5) defined, one simply choose y = C = (1 0)e (or y = ei) as 
the output of the error system. In the synchronization approach, one writes 
y = h(e) = e\ = X\m — X\s, and consequently the error system (5) can be 
written in the general form 



X = f{e)+g{e)U . 



(7) 



The error system (J2J) should be stabilyzed at the origin or in an arbitrary 
small neighbourhood of it. More details on the synchronization conditions 
are provided in the papers |3|4| that we employ to obtain the control function 




U 



-1 

7 



(P + 8) , 



(8) 



where 7 and S are real- valued functions obtained by means of Lie derivatives 
of h{e) as follows 



7 



L g L p f - l h(e) 



P = L p f h(e) 



(9) 



where p is a positive integer that determines the so-called relative degree of 
the system (see On the other hand, the desired dynamics, i.e., directed 
towards the origin, is dictated by 



5 = Kp\e\ + Kp 2 e2 + Kp^e^ . 



(10) 



Thus, performing the Lie derivatives and regrouping the terms, we obtain 
the function U of the form 



U 



-1 

KT 



(-K 2 e'i + KTe' 3 ) + Kp\e\ + Kp 2 e 2 + K P3 e 3 



Using the change of variables = JQjvf — Xips, where the latter vector is 
the column vector formed by the triad T, N, and B, (i = 1, 2, 3), the control 
can be written as a function of the states of the two systems 



U = 



-1 

KT 



■ k 2 (X 2 - kN) - kt (-aX 3 - X x + X X X 2 ) + (tN) 



+ 5 



where 5 = K P1 {X 1 - T) + K P2 (X 2 - N) + K P3 (X 3 



(12) 

B). Notice that 



7 = LgLjrh(x) = kt is a nonzero constant. Therefore the control signal is 
defined for any T,N,B, X\,X 2 , and X 3 . In addition, one should choose the 
values of the constants Kpi in such a way that the differences in 5 go to 
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zero. Applying the dynamics generated by ()12|) leads to the synchronization 
matrix 



From this synchronization matrix one can see that the first two states of 
both systems are synchronized. However, the state B = X%m + kXim is 
extended by the term kX\m, i.e., the state B is the sum of two states of the 
Malasoma oscillator; since the latter is chaotic, one concludes that the state 
B is also chaotic. 

We display the phase locking between the corresponding phases of the 
two oscillators in Figs. 3 and 4 where the phase locking of the states T and 
X\ and N and X2, respectively, shows that the two pairs of states are syn- 
chronized. In Fig. 5, we see that the B and X3 states are not synchronized. 
Thus, following the terminology of [H], we are in the situation of a general- 
ized synchronization. In Fig. 6, the two already synchronized systems are 
shown in the three dimensional space. One can notice that the FS system is 
'above' the Malasoma oscillator, and that the two systems are in a chaotic 
phase. Finally, in Fig. 7, the control signal used to achieve the generalized 
synchronization of this paper is displayed. 

In summary, we have shown here in a concrete way how the simple chaotic 
dynamics of Malasoma type can be imposed to the linear Frenet-Serret evo- 
lution of space curves. 
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Figure 1: Three dimensional dynamics of the Frenet-Serret system with 
k = 1 and r = 0.9. 
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Figure 2: Three dimensional dynamics of the Malasoma system with a = 
2.025. 
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Figure 3: Phase locking of the T states of the FS system to the X\ states of 
the Malasoma system. 
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Figure 4: Phase locking of the N states of the FS system to the X2 states 
of the Malasoma system. 
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B 

Figure 5: Phase behaviour of the states B and X3 of the FS system and 
Malasoma system, respectively. 
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Figure 6: Plot of the two chaotic attractors once they are synchronized in 
the generalized form. 
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Figure 7: The control signal U imposed to the Frenet-Serret system in order 
to display Malasoma chaos. 



